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1.  INTRODUCTION 


Given  an  n  by  n  matrix  A,  n-vector  b  and  a  mapping  m  from  Rn  into  Rn, 
the  implicit  complementarity  problem  (ICP) ,  denoted  by  the  triple  (A,b,m) 
is  to  find  a  vector  x  such  that  the  conditions  below  are  satisfied: 

(1)  Ax  +  b  >  0,  x  >  m(x)  and  (Ax  +  b)^(x  -  m(x))  =  0. 

If  m  is  the  constant  mapping  m(x)  =  c,  then  the  problem  (1)  becomes 

Ax  +  b  >  0,  x  >  c  and  (Ax  +  b)T(x  -  c)  =  0 
which  clearly  is  equivalent  to  the  linear  complementarity  problem  (LCP) 

Ay  +  (b  +  Ac)  >0,  y  >  0  and  yT  [Ay  +  (b  +  Ac)  ]  =  0 

under  the  translation  y  =  x  -  c.  If  the  matrix  A  is  nonsingular,  then  the 
ICP(l)  is  equivalent  under  the  identification  y  =  Ax  +  b,  to  the  nonlinear 
complementarity  problem  (NLCP) 

(2)  y  >  0,  g(y)  >  0  and  yTg(y)  =  0 

with  g(y)  =  -A  *b  +  A  1y  -  m(-A  H  +  A  ^y) .  Conversely,  it  is  obvious 
that  an  arbitrary  NLCP (2)  can  be  written  in  the  form  (1)  with  A  =  I,  b  =  0 
and  m(x)  *  x  -  g(x). 

More  generally,  if  the  mapping  m  is  piecewise  linear,  i.e.  if 

m(x)  =  max  [A*x  -  b^] 

1  <  i  <  k 

where  each  A^  and  each  b^  are  n  by  n  matrix  and  n-vector  respectively,  then 


the  ICP(l)  becomes  the  generalized  linear  complementarity  problem  (see 
Cottle  and  Dantzig  [9  ]  and  Mangasarian  [21]) 

Ax  +  b  >  0,  (I  -  A1)  x  +  b1  >  0  for  each  i  -  1 . k 

k 

TT  i  i 

(Ax  +  b)^  ^  [(I  -  A)  x  +  b  ]j  =«0  for  each  j  =  l,...,n 

which  obviously  is  equivalent  to  the  problem  of  finding  a  solution  to 
the  system  of  piecewise  linear  equations  (see  Eaves  [14]  e.g.) 

min  {Ax  +  b,  min  [(I  -  aS  x  +  b1-]}  =  0. 

1  <  i  <  k 

The  ICP(l)  has  been  studied  in  Dolcetta  [13]  as  a  discretized  version 
of  a  continuous  variational  problem.  There,  a  certain  implicit  (point-to- 
set)  mapping  was  defined  and  used  to  analyze  the  problem. 

The  present  paper  is  the  first  of  a  two-part  research  on  the  ICP. 

Here,  we  study  the  (global)  existence  and  uniqueness  of  solution  to  the 
ICP  and  investigate  the  connection  of  the  ICP  with  a  least-element  theory. 

In  the  second  part  of  this  research  (which  is  presented  in  an  accompanying 
paper  [28]),  we  establish  a  convergence  theory  for  a  certain  iterative 
algorithm  to  solve  the  ICP  and  show  how  this  theory  unifies  and  extends  many 
existing  convergence  results  for  the  successive  overrelaxation  methods  to 
solve  a  LCP.  The  fundamental  device  employed  throughout  is  a  certain 
implicitly  defined  mapping  F  which  is  constructed  from  a  splitting  of  the 
matrix  A.  The  construction  of  F  is  largely  based  on  ideas  extending  those 
used  in  Dolcetta 's  work  [13]  and  Aganagic's  recent  doctoral  dissertation  [1], 
-  The  organization  of  the  remainder  of  the  paper  is  as  follows.  Section 
2  contains  two  parts.  In  the  first  part,  we  shall  establish  some  preliminary 


results  about  the  I CP.  In  particular  we  shall  derive  a  basic  relationship 
between  the  ICP  and  a  certain  variational  problem.  In  the  second  part, 
we  shall  introduce  the  mapping  F,  lay  down  a  fundamental  relationship 
between  F  and  solution  to  the  ICP  and  review  some  basic  facts  which  are 
useful  in  the  subsequent  study.  Section  3  is  concerned  with  the  existence 
and  uniqueness  of  solution  to  ICP.  Our  main  goal  is  to  establish  some 

f 

sufficient  conditions  (on  the  matrix  A  and  the  mapping  m)  under  which  the 
ICP(A,b,m)  will  have  a  unique  solution  for  all  vectors  b.  The  conditions 
derived  are  closely  related  to  those  obtained  in  Megiddo  and  Kojima  [22] 
for  the  NLCP.  The  approach  taken  is  based  on  ideas  used  in  Aganagic  [1], 
Eaves  [14]  and  Megiddo  and  Kojima  [22].  In  Section  4,  we  shall  study  the  ICP 
in  connection  with  a  least-element  theory.  In  particular,  we  present  some 
conditions  under  which  the  ICP  will  have  a  solution  characterized  as  the 
least-vector  of  the  feasible  set,  and  establish  the  monotone  convergence  of 
a  certain  iterative  scheme  (defined  by  the  mapping  F)  to  solve  the  ICP. 

Some  results  of  this  last  section  have  previously  appeared  in  the  author's 
dissertation  [25]. 
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2.  PRELIMINARY  DISCUSSION 

2.1.  Some  Basic  Facts.  It  is  well-known  that  there  is  a  fundamental 
connection  between  the  NLCP(2)  and  the  variational  inequality  problem. 
The  latter  is  to  find,  given  a  set  K  in  Rn  and  a  mapping  g  from  R° 
into  itself,  a  vector  u€K  such  that 

(x  -  u)T  g(u)  >  0  for  all  x€k. 

See  Cottle  [7]  and  Karamardian  f 16 ] .  In  the  proposition  below,  we 
establish  an  extension  of  this  connection  for  the  ICP. 

Proposition  2.1.  A  vector  u  solves  the  ICP(l)  if  and  only  if  it 
satisfies  the  conditions 

(3)  u€K(u)  and  (x  -  u)T  g(u)  >0  for  all  x€K(u) 
where  g(x)  is  the  affine  mapping 
g(x)  -  Ax  +  b 
and  K(u)  is  the  set 

K(u)  *  Cx€Rn  :  x  >  m(u) } . 

The  problem  of  finding  a  vector  u  satisfying  (3)  is  known  as  a 
quasi -variational  inequality  problem  (QVI) .  This  terminology  was 
due  to  A.  Bensoussan  and  J.  L.  Lions.  In  fact  in  a  series  of  papers 
(2,  3,  4,  5]  these  authors  have  defined  both  the  QVI  and  the  ICP  in  an 
infinite-dimensional-space  setting  and  showed  how  various  impulsional 
control  and  optimal  stopping  problems  can  be  formulated  as  (infinite- 


dimensional)  QVI's  and  ICP's.  Here  we  shall  not  go  into  further  dis¬ 
cussion  of  this  matter,  but  simply  refer  the  interested  readers  to  the 
cited  references  for  more  details.  We  would  like  to  point  out  however, 
that  these  impulsional  control  and  optimal  stopping  problems  in  their 
discretized  versions,  have  provided  another  source  of  application  of 
the  ICP  defined  in  (1).  See  Dolcetta  [13], 

Proof  of  Proposition  2.1.  Suppose  that  u  satisfies  (3).  Then  u  >  m(u) 
T 

and  (m(u)  -  u)  g(u)  >  0  as  m(u)  €K(u) .  Moreover,  2u  -  m(u)  €k(u)  so 
that 

0  <  (2u  -  m(u)  -  u)T  g (u)  =  (u  -  m(u))T  g(u) . 

X  n 

Hence  (u  -  m(u))  g(u)  =  0.  Finally,  if  x€r+,  then  x  +  u€k(u) .  Thus, 

0  <  (x  +  u  -  u)T  g (u)  =  xTg(u). 

Consequently,  g(u)  >  0.  This  shows  that  u  solves  the  ICP(l). 

Conversely,  suppose  that  u  solves  the  ICP(l).  Clearly  u€k(u)  .  If 
x€k(u),  then  x  >  m(u) .  Hence,  we  have 

(x  -  u)T  g (u)  =  (x  -  m(u))Tg(u)  +  (m(u)  -  u)T  g(u) 

*  (x  -  m<u))Tg(u)  >  0. 

In  other  words,  u  satisfies  (3).  Q.E.D. 

An  Immediate  consequence  of  the  proposition  is  that  one  may 
formulate  the  ICP  as  a  fixed  point  problem.  More  specifically,  we  have 
Corollary  2.2.  A  vector  u  solves  the  ICP(l)  if  and  only  if  u€k(u)  where 

K(u)  ”  [v^K(u)  :  v^(Au  +  b)  *  min  x^(Au  +  b)}. 

x€K(u) 
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It  is  a  known  fact  that  the  NLCP(2)  has  at  most  one  solution  if  g 

1  2 

is  a  P-function.  i.e.,  if  for  every  vectors  y  ^  y  ,  it  holds  that 

max  (y1  -  y2)  (gCy1)  -  g(y2)).  >  0. 

1  <  i  ^  n  1 

See  More  [23]  e.g.  Specializing  this  to  the  ICP(l),  we  obtain 

Proposition  2.3.  Suppose  that  the  mapping  m  and  the  matrix  A  satisfy 

1  2 

the  condition:  for  every  x  /  x  , 

(4)  max  [(x1  -  mCx1))  -  (x2  -  m(x2))  ]  (A(x1  -  x2))  >  0. 

1  <  i  <  n 

S  2= 

Then  the  ICP(l)  has  at  most  one  solution  for  every  vector  b. 

Proof.  In  fact,  the  given  condition  (4)  implies  that  the  matrix  A  is 
nonsingular.  Hence,  as  pointed  out  in  the  Introduction,  the  ICP(l)  is 
then  equivalent  to  the  NLCP(2)  where  the  mapping  g  is  given  by 

g(y)  =  a"1^  -  b)  -  m(A  1(y  -  b))  . 

It  is  easy  to  deduce  that  condition  (4)  is  equivalent  to  the  fact  that 
g  is  a  P-function.  Thus  the  desired  conclusion  follows  readily.  Q.E.D. 

Note  that  the  condition  (4)  is  not  sufficient  to  imply  the 
existence  of  solution  to  the  ICP. 

2.2.  Definition  of  F.  A  splitting  of  the  matrix  A  is  a  representation 

(5)  A  -  B  -  C. 

A  P-splitting  (Z-splitting)  is  one  where  B  is  a  P-matrix  (Z -matrix) . 

Recall  that  a  P-matrix  (Z-matrj.x)  is  a  real  square  matrix  whose 


principal  minors  are  all  positive  (whose  off-diagonal  entries  are  all 
non-positive).  A  K-splitting  is  one  where  B  is  a  K-matrix  (i.e.  both 
a  P-  and  Z -matrix).  The  trivial  splitting  is  where  B  is  the  identity 
matrix.  Note  that  the  trivial  splitting  is  a  K-splitting.  We  shall 
denote  the  splitting  (5)  by  the  pair  (B,C). 

Given  a  P-splitting  (B,C)  of  the  matrix  A,  we  define  a  mapping  F 
from  Rn  into  Rn  as  follows.  For  each  vector  u€Ra,  F(u)  is  the  unique 
solution  to  the  LCP 

Bx  +  (b  -  Cu)  >0,  x  >  m(u)  and  (x  -  m(u))T(Bx  +  (b  -  Cu))  =0. 

Since  B  is  a  P-matrix,  F(u)  is  well-defined  (see  Samelson  et  al  [30]). 

In  general,  the  mapping  F  depends  on  the  splitting  and  possesses  no 
explicit  representation  (except  in  certain  special  cases).  However, 
each  F(u)  can  be  computed  by  a  number  of  efficient  algorithms,  see 
Cottle  [8],  If  (B,C)  is  the  trivial  splitting,  then 
F(u)  =  min(m(u),  (I  -  A)u  -  b) . 

We  point  out  that  if  m  is  the  zero  mapping,  then  the  latter  F  is  precisely 
the  mapping  used  in  Aganagic  [1]  to  study  the  LCP.  If  (B,C)  is  a  K- 
splitting,  then  F(u)  has  the  characterization  below.  Moreover,  it  can  be 
computed  by  a  special-purpose  algorithm  due  to  Chandrasekaran  [6], 

Theorem  2.4.  (Cottle  and  Veinott  [12]).  Let  (B,C)  be  a  K-splitting.  Then 
F(u)  is  characterized  as  the  least-element  of  the  polyhedral  set 

X (u)  «  [x  6  Rn  :  Bx  +  (b  -  cu)  >0,  x  >  m(u)}. 

Moreover,  it  is  the  unique  solution  to  the  linear  program 
T 

min  e  x  subject  to  x  €  X(u) 


where  e  is  any  positive  vector 
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Recall  that  the  least-element  x  of  a  set  S  in  R°  is  one  such  that 
x  <  x  for  all  x€s.  Another  P-splitting  that  is  of  particular  interest  is 
the  one  where  B  is  a  lower  or  upper  triangular  matrix  with  positive 
diagonal  entries.  In  fact,  it  is  easy  to  verify  that  in  this  case,  the 
components  of  F(u)  are  given  inductively  by  :  for  each  i, 

-  F(u>  -  max  [m(u)  ,  -  [(b  -  Cu)  -  Z  B  F(u)  J/B  } 

i  i  ij<iijjn 

or 

F(u)  =  max  fm(u)  ,  -  [(b  -  Cu)  -  Z  B  F(u)  ]/B  } 

depending  on  whether  F  is  lower  or  upper  triangular. 

If  the  mapping  m  is  piecewise  linear  (as  in  the  case  of  the  gener¬ 
alized  linear  complementarity  problem),  then  so  is  the  mapping  F.  The 
proposition  below  asserts  a  fundamental  relationship  between  solution(s) 
to  the  ICP(l)  and  the  mapping  F.  Its  proof  is  easy  and  is  thus  omitted. 

Proposition  2.5.  A  vector  u  is  a  solution  to  the  ICP(l)  if  and  only 
if  it  is  a  fixed  point  of  the  mapping  F  for  any  P-splitting  of  the 
matrix  A.  In  particular,  if  some  F  has  a  unique  fixed  point,  then  the 
ICP(l)  has  a  unique  solution. 

We  state  an  important  property  for  a  K-matrix  in  the  theorem 
below.  A  proof  can  be  found  in  the  reference  cited. 

Theorem  2.6.  (Fiedler  and  Ptak  [15]).  If  B  is  a  K-matrix,  then  so  is 
each  of  its  principal  submatrices.  The  inverse  B  *  exists  and  is  non¬ 
negative. 

We  close  this  section  by  noting  that  the  (point-to-set)  mapping 
used'  in  Dol'cetta  [131  to  study  the  ICP  corresponds  to  the  mapping  F 
defined  by  the  splitting  (B,0)  with  B(*A)  being  a  Z-matrix.  More  about 
this  will  be  said  in  Section  4. 


3.  EXISTENCE  AND  UNIQUENESS 


Associated  with  the  mapping  F  defined  by  a  P-spLitting  of  the  matrix 
A,  we  define  the  mapping 

G(x)  -  x  -  F(x) . 

In  view  of  Proposition  2.5,  we  have 

Theorem  3.1.  A  sufficient  condition  for  the  ICP(l)  to  have  a  unique 
solution  is  that  the  mapping  G  is  a  bijection  for  some  P-splitting  of 
the  matrix  A. 

Our  goal  in  this  section  is  to  derive  some  sufficient  conditions 
(on  the  function  m  and  matrix  A)  so  that  for  each  vector  b,  the  mapping 
G  is  a  global  homeomorphism  and  is  thus  a  bijection.  This  then  implies 
that  the  ICP(l)  will  have  a  unique  solution  for  each  vector  b. 

This  approach  of  using  a  global  homeomorphic  mapping  to  study  the 
(global)  existence  and  uniqueness  of  solution  to  a  complementarity 
problem  has  previously  been  taken  by  Megiddo  and  Kojima  [22]  in  the 
case  of  the  NLCP(2).  It  is  not  clear  to  the  author,  however,  whether 
the  specific  mapping  g(y+)  +  y”  where  y+  and  y  are  the  positive  and 
negative  part  of  the  vector  y  used  in  the  reference  can  be  obtained  by 
our  approach. 

The  rest  of  this  section  is  divided  into  two  subsections.  In 
the  next  subsection,  we  shall  first  derive  conditions  under  which  the 
mapping  G  is  a  local  homeomorphism.  In  Subsection  3.2  we  shall  study 
the  global  property  of  G. 


10 


3.1.  Local  Properties.  Throughout  the  rest  of  this  section,  we  let 
the  mappings  F  and  G  be  defined  with  respect  to  a  fixed  but  arbitrary 
K-splitting  (B, C)  of  the  matrix  A. 

Proposition  3.2.  If  the  mapping  m  is  convex,  then  so  is  F. 

1  2 

Proof.  For  any  two  vectors  u  and  u  and  scalar  X€[0,1],  the  vector 

XF(u^)  +  (l-X)F(u^)  belongs  to  the  set  X(u)  where  u  =  Xu*'  +  (1  -  X)u^. 

As  F(u)  is  the  least  element  of  this  latter  set  (Theorem  2.4),  it 
-1  2 

follows  that  F(u)  <  XF(u  )  +  (1  -  X)F(u  ),  proving  the  convexity  of 
F.  Q.E.D. 

Recall  that  if  f  is  a  convex  mapping  from  R°  into  Rn,  then  an  n 
by  n  matrix  E  is  a  subgradient  of  f  at  the  point  x  if 

f(y)  >  f00  +  E(y  -  x)  for  all  y€Rn. 

In  particular,  each  row  E^  of  E  is  a  subgradient  of  the  i-th  com¬ 
ponent  function  f^  at  the  point  x.  We  denote  by  3f(x)  the  set  of 
subgradients  of  f  at  x.  If  a  is  a  subset  of  {l,...,n},  then  by  3fa(x) 
we  denote  those  rows  of  3f(x)  indexed  by  a.  In  particular,  each  af^x) 
is  a  set  of  row  vectors. 

Returning  to  our  mapping  F,  we  define  for  each  vector  u,  three 
partitioning  subsets  of  [l,...,n]: 

<*,  00  ■  (i,  :  F(u),  =  m(u) 

l  1  h  Li 

a  (u)  -  [i  :  F(u)  >  m(u) 

L  2  2 

a,  00  *  £i,  :  F(u)  =*  m(u) 

j  i  i3  13 

For  each  subset  v(u)  of  a^OO*  we  let  £(u)  be  its  complement  in  a^(u) 


and  (BF(u)  +  (b  -  Cu))  >0} 

1 

and  (BF(u)  +  (b  -  Cu))  *  0} 

12 

and  (BF(u)  +  (b  -  Cu))  =  0}  . 

*3 
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and  let  Y(u)  *  ar^u)  U  Y(u).  If  the  point  of  reference  u  is  obvious  in 
the  context,  we  shall  omit  it  from  these  index  sets.  For  each  subset 
Y  of  {l,...,n},  we  shall  denote  (B  )  *  by  B^.  By  Theorem  2.6,  each 
B^  is  nonnegative. 

Proposition  3.3.  Let  m  be  convex  and  let  u€rD.  Then 


(6a)  SF  (u)  =  3m  (u) 
1  1 


(6b)  (u)  3  [-  _S  AS  -  [B- 


5m 


Y  c  «,(u)  V  «2Y  Y«^  UP^U  P 


,(u)  -  C- ]  : 


A  6  A  >  0  and 

•  K!  = f 


Y  C  or2(u)  Y 


1} 


(6c)  For  each  i^  €  Qf^(u), 


i3(u)  =  i3e  Y  c«3<u)  Vi3Y  [BV«  U^au  P(u)  “  ty 

+  (1  -  X^)Smi  (u)  ;  A  ,  \ ^  €  R;  \  >  0  and 


2  A  +  Xy  m  l] 

i  6  Y  f  af3(u)  Y 

where  in  (6b)  and  (6c),  C-  denotes  the  rows  of  the  matrix  C  indexed  by  y. 
To  prove  the  proposition,  we  need  the  following  lemma. 


Lemma  3.4.  The  vector  h  is  a  subgradient  of  a  convex  function  f  :  R°  -*  R 

at  the  point  x  if  and  only  if  there  is  a  neighborhood  Nx  of  x  such  that  for 

all  y  €  N 
J  x 

(7)  f (y)  >  f(x)  +  hT(y  -  x). 

The  lemma  says  that  the  vector  h  is  a  subgradient  if  and  only  if  the 
subgradient  inequality  (7)  is  satisfied  locally. 


^ 
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Proof  of  Lemma  3.4.  Only  the  sufficiency  part  needs  to  be  proved.  Let  z 
be  an  arbitrary  point  in  Rn.  Then  there  is  a  X  6  (0,1)  such  that  the 
vector  y  ■  \z  +  (1  -  X)x  belongs  to  the  neighborhood  N^.  By  the  monotonicy 
property  of  a  convex  function  (see  Stoer  and  Witzgall  [31,  p.  135])  we  have 

f(z)  -  f(x)  >  (f(y)  -  f(x))/X  >  hT(y  -  x)/X  =  hT(z  -  x) 


where  the  second  inequality  follows  from  the  assumption.  Q.E.D. 


Proof  of  Proposition  3.3.  Since  both  m  and  F  are  convex,  they  are 
continuous  (see  Stoer  and  Witzgall  [31,  p.  135]).  Hence,  there  exists  a 
neighborhood  Nu  of  u  such  that  for  all  v  €  N^,  we  have 

[BF(v)  +  (b  -  Cv)]a  >0  and  (F(v)  >  m(v))a  . 

By  complementarity,  it  follows  that  for  all  v  6  Nu> 


(F  (v)  *  m(v)) 


«x(u)  ' 


Consequently,  (6a)  follows  from  the  last  lemma.  To  prove  (6b) ,  observe 
that  for  each  y  c  (u) ,  we  have 

<F<u»;  -  y  B  H u  t  *  0.  -  CU)-I 

and  for  all  v 

[BF (v)  +  (b  -  Cv) ]-  >  0. 

Thus,  it  follows  that. 


1  -  \=>  U  B  F(v,«1  U  B  +  ^ 

2  ‘  V  U  B  ”(V>«1  u  B  +  <b  '  °,)?1 

S  •  B«  y  tBv«  u  B  <m(u)  +  E(v  ■  “^U  B  +  <b  *  ^Y1 


-13 


where  E  .  a(u) .  Note  that  we  have  used  the  nonnegativity  of 

Up  a^  U  P 

the  matrix  -B  -B-  . .  a  to  get  the  last  two  inequalties.  Hence,  for 

a2Y  YQfl  U  * 

each  vector  v 


F<v>„2  >  F(»>^  -  u  p  *Bi  u  ,  -  Oy]  <*  -  »>• 

Assertion  (6b)  now  follows  from  the  convexity  of  the  subgradient  set 

3F  (u) .  The  last  assertion  (6c)  can  be  proved  in  a  similar  fashion. 
a2 

Q.E.D. 

Corollary  3.4.  Let  m  be  convex  and  u  6  Rn .  Then 


(8a) 

(8b) 


(8c) 


^ 

Vu)  -  \ +  v  1 .3<o)  MvS 

j.y  e  Rl*2l*la4l,  Ay  >  0 


For  each  i^  6  or^u) 


U  8  am«1  U  3(u)  "  Cv  : 

and  S'  A  =  i} 

V  <=  «3(u)  Y 


3GS(U)  2  \  +  i  s  v  =  «3w  'S.  Uf>«iu  s<u)  '  ^ 

-  (1  -  \^)am  (u)  :  \  ,  6  R  ;  \  ,  \  >  0 

md  i  €  Y*<=  «3(u)  \  +  "  1}* 


Here  I  denotes  the  rows  of  the  identity  matrix  indexed  by  y. 

A  remark  is  in  order.  If  m  is  convex,  then  the  mapping  G  is  con¬ 
cave.  A  subgradient  of  a  concave  function  is  a  vector  whose  negative 
is  a  subgradient  of  the  negative  of  the  function  (which  is  convex). 

It  has  been  shown  in  Eaves  [14 ]  and  Aganagic  [1 ]  that  if  the 
mapping  m  is  piecewise  linear,  then  the  mapping  G  defined  with  respect 


mm 
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to  the  trivial  splitting  of  A  is  locally  univalent  at  u  provided  that 
each  subgradient  matrix  in  3G(u)  is  nonsingular.  The  proposition  below 


extends  this  result. 


Proposition  3.5.  Let  m  be  convex  and  u  €  R  .  Suppose  that 

(9)  For  each  E  €  3m (u) ,  and  each  (A^,A^)€  r  I  a2 ! x !  a2 1  x  rI^M0^! 
satisfying 

(10)  0  <  A  .  E  A  =  I  and  0  <  A'  E  A'  <  I, 

~  Y  Y  =  a3  Y  *  Y  ’  V  <=  or3  (u)  Y  " 


the  matrix  below  is  nonsingular 


H  ■  \ 


\  +  yi  n3(u)  Vv  [BY«1  CY] 


1  +  z 

*3  YC«3(u) 


\Baf3Y  ^YOfj  US^UrV*  (I 


Z  A*)E 

Y  c  a3(u)  Y  a3 


Then  the  mapping  G  is  locally  univalent  at  u. 

Proof.  By  [Thm.  24.7,  29],  each  dm^u)  is  a  compact  set.  Thus  so  is  dm(u)  ■ 

Sm- (u)  x  ....  x  dm  (u) .  Hence,  there  exists  an  e  >  0  such  that  if  N  is  the 
1  n 

unit  ball  in  the  space  of  n  by  n  matrices  (considered  as  the  product  space  of 
n  copies  of  R  ),  then  the  matrix  E  is  nonsingular  for  every  E  c  H  +  *N  for 
some  matrix  H  of  the  form  (11). 


Suppose  that  G  is  not  locally  univalent  at  u.  Then  there  exist  sequences 
[x^]  and  Cy1}  both  converging  to  u  such  that  for  each  i,  G(xi)  =  G(yS  and 
tl~  ^  .  With  no  loss  of  generality,  we  may  assume  that  the  two  sequences  are 

chosen  so  that  for  each  i. 


15 


[BFCx1)  +  (b  -  Cx1)  ]  (u)  >  0,  [BF(yi)  +  (b  -  Cy1)  ]  j  >  0 


(F(xS  >  mCx1)) 


<*2(u) 


,  (F(,S  >»(y1))<l2(u). 


Let  5  >  0.  By  [Thm.  24.5,  29],  there  is  an  index  j  such  that  dm(x^)  U  Bm(y^) 
C  3m(u)  +  *N.  Let  E^  €  Sm(x^)  and  E^  €  Sm(y^) .  Define  four  mutually 
disjoint  subsets  of  ar3(u)  as  follows: 


^(u)  =  Cix  €  0f3(u) 

:  F  (xj)  =  m  (xj) 

1  \ 

and 

F.  (y-5)  =  m  (yj)} 
X1  L1 

P2(u)  =  [i2  €  c*3(u) 

:  F  (xJ)  =  m  (xj) 

2  2 

and 

F.  (y-5)  >  m.  (yJ)  ] 

2  i2 

P3(u)  =  [i3  €  ar3(u) 

:  F.  (x^)  >  m.  (x^) 

L3  x3 

and 

f.  (yj)  =  m  (yj)} 
3  3 

-  3 

P4(u)  =  a3(u)  \ 

^(u) . 

Now  for  each  i  €  P^(u)  U  Qf^(u),  we  have 

x|  -  mi(xj)  -  Gi(xJ)  =*  Gi(yJ)  -  yj  -  m^) 

By  the  fact  that  (E*,  E^)  €  5m^(x^)  x  Sm^Cy^),  it  follows  that 


(It  -  E*)(yj  -  xj)  >  0  and  (I.  -  E‘)(xJ  -  yJ) 


>  0 


1  2 

where  1^  is  the  i-th  row  of  the  identity  matrix  and  where  E^  and  E^  are 

1  2 

the  i-th  rows  of  the  matrices  E  and  E  respectively.  Hence,  there 
exists  a  X.^  €[0,1]  such  that  for  E^  =  (1^  -  E^)  *  X^CI^  “  eJ)  +  (1  •  X^) 


(Ii  -  E^),  we  have 


(12)  Et(xj  -  yJ)  -  0. 


Note  that  Ei  €  3m^(u)  +  SN.  Moreover,  for  i  €  B2(u),  we  have 
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-  m^x3)  -  Gt(xJ)  =  G1(yj)  <  y|  -  m^y-3). 

Thus, 

(13)  (It  -  E^)  (yj  -  xj)  >  0. 


On  the  other  hand,  we  have  by  complementarity, 

[BF(yj)  +  (b  -  Cyj)]i  -  [B(y3  -  G(y3))  +  (b  -  Cy3)^  =  0 

<  [BF(x3)  +  (b  -  Cy:3)]i  =  [B(x3  -  G(x3))  +  (b  -  Cx3)  ].^ 
Thus,  it  follows  that 

(C±  -  Bi)  (xj  -  y3)  <  0. 


Similarly,  we  can  show  that  for  each  i  €  ^(u)» 


(Ii  “  EP  >  0  and  (ci  ■  Bi>  (yJ  *  *J)  i  °- 

Moreover,  for  each  i  €  (J^(u)  U  we  ^ave 

(Ci  -  B^  (x3  -  y3)  -  0. 

Now,  for  an  index  i  €  *-et  Vj_  *  13  I^OO*  T^en 

**1 + 5iv1v1 » h\ « \ + \*A  -  v  <iJ  • yJ) 

i  &i  +  5ty1tBv, 

-  [Ii  +  ®iyt  t\> 1  (xJ  '  yJ)  ‘  °’ 

Since  (Ii  -  E^)  (x3  -  y*-)  <  0  by  (13),  there  exists  a  €[0,1]  such  that 
for  Et  -  It  +  U  &1Ecr1  U^+  B?iP3EP3  '  CV1I  "  (1  "  Xi)Ei* 

(12)  holds. 


J  -  J, 


’iVVl  U  ^  U  P,  +  \*3h2  ’  BYt]1  (XJ  -  yJ> 


Similarly  for  i  €  6^(0)  let  *  0^(u)  ^  B^(u).  Then  there  exists  a 


Xi  €(0'U  s"'h  th“  for  El  ‘  rt  +  Vlv^Oj  U  Bj'oj  UB/  8vtB2EB2  ‘  \] 

2 

-  (1  -  X^E^,  (*2)  holds  as  well.  Finally,  for  i  €  3^(u)  Uar2(u), 
let  *  (^(u)  U  311(1  =  U  &.j(u).  Then  there  exists  a 

xi S[<M1  such  tha,:  £or  Ei  *  ri +  1i5iy1'v1  u  \\ u  h +  *Wh '  % 

*  (1  -  V5i«1  u  B1Ec1  U  s2  +  \B2EB2  •  %'•  (I2)  holds- 


Let  E  be  the  matrix  whose  rows  are  as  defined  above.  By  choosing 

6  >  0  small  enough,  we  have  E  €  H  +  cN  where  H  is  a  matrix  of  the  form  (11) . 

~  j  ]  ~ 

Since  E(x  -  y  )  =  0,  we  obtain  a  contradiction  to  the  nonsingularity  of  E. 

Q.E.D. 


Remark.  By  Corollary  3.4,  each  matrix  given  in  (11)  is  a  subgradient  in 
dG(u) . 

A  mapping  f  from  Rn  into  itself  is  said  to  be  inverse  locally 
Lipschitz  continuous  at  the  point  x  if  there  is  a  neighborhood  N-x  of  x 
and  a  constant  scalar  k  such  that 


|x  -  yfl  <  k  |ff(x)  -  f(y)||  for  all  y  €  N^. 


It  was  shown  in  Aganagic  [1]  that  in  the  case  of  the  LCP,  the  mapping 
6  defined  with  respect  to  the  trivial  splitting  of  the  matrix  A  is 
always  inverse  locally  Lipshitz  continuous  provided  that  certain 


-  'HRs'  - 
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regularity  conditions  are  satisfied.  In  what  follows,  we  extend  this  re¬ 
sult.  We  define  for  each  subset  y  of  {l,...,n},  a  mapping  C*  from  Rn  into 
itself  by 

(•1Y(x))y  -  A^x  and  Ot«Y(x))p  “  x^  -  m^(x) 
with  8  the  complement  of  y. 

Proposition  3.6.  Suppose  that  for  each  subset  y  of  {l,...,n},  the  mapping 

y 

L  defined  above  is  inverse  locally  Lipschitz  continuous  at  the  point  u. 
Then  so  is  G. 

Proof.  Let  Ny  be  a  neighborhood  of  the  point  u  so  that  for  each  y  €  N^, 
the  following  holds 

[BF(y)  +  (b  -  Cy)]  (u)  >  0,  <F(y)  >  m(y))^(u) 

and  for  all  subsets  y  containing  ar^(u) , 

lfu  *  yll  <  kH‘tV(«)  -  ^Y(y)II 

for  some  scalar  k.  Let  y  be  an  arbitrary  point  in  and  define  two 
index  sets  8'  and  y'  by 

8'  *  {i  €  <*3(u)  ;  F(y)i  «  m(y)i},  y'  »  fi  €  «3(u)  :  F(y)£  >  m(y)t] 

By  letting  y  ■  «2(u)  U  y*  and  8  =  «^(u)  U  8',  we  have 
[G(u)  -  G(y)]p  *  (u  -  m(u))p  -  (y  -  m(y))^ 

[B(u  -  G(u) )  +  (b  -  Cu) ]  -  [B(y  -  G(y))  +  (b  -  Cy) 1  -  0 

or  equivalently, 

ly(u)  -  *Y(y)  -  /  Ip  \ 


(G(u)  -  G(y)) 


Consequently,  for  each  y  €  N^, 

l|u  -  y||  <  K||G(u)  -  G(y) || 

where 

) 11 . 

Therefore,  the  mapping  G  is  inverse  locally  Lipschitz  continuous  at  the 
point  u.  Q.E.D. 

It  is  easy  to  see  that  if  the  mapping  m  is  a  constant,  a  sufficient 

y 

condition  for  each  mapping  L  to  be  inverse  locally  Lipschitz  continuous 
(everywhere)  is  that  the  matrix  A  is  nondegenerate,  i.e.  each  principal 
submatrix  of  A  is  nonsingular.  This  nondegeneracy  of  A  is  precisely 
the  regularity  property  used  in  Aganagic  [1]  to  establish  Proposition  3.6 
for  the  LCP  in  the  case  of  the  trivial  splitting.  We  remark  that  no 
convexity  assumption  on  the  mapping  m  is  needed  in  order  for  this  last 
proposition  to  hold. 

The  next  theorem  gives  sufficient  conditions  for  G  to  be  a  local 
homeomo  r ph i  sm . 

Theorem  3.7.  Let  m  be  convex  and  u  €  Rn.  Suppose  that  condition  (9) 
is  satisfied  at  u.  Then  the  mapping  G  is  a  local  homeomorphism  mapping  a 
neighborhood  of  u  onto  a  neighborhood  of  G(u) . 

Proof.  Proposition  3.5  implies  that  there  is  a  neighborhood  Ny  of  u  where  G 
is  injective.  Since  G  is  continuous,  by  [Theorem  A.l,  22],  G  maps  N 
homeomorphically  onto  the  neighborhood  G(Nu>.  Q.E.D. 


0  ?  or. 


max 
<u),  Y 


of2(u) 


Remark.  The  inverse  continuity  of  G  is  a  consequence  rather  than  an 
assumption  of  Theorem  3.7. 
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3.2.  Global  Homeomorphism.  To  continue  the  discussion,  we  say  that  a 
mapping  f  :  Rn  -*  Rn  is  norm-coercive  if  for  every  sequence  [xk]  with 
-•  «,  we  have  [|f(xlc)|[  -*  ». 

Proposition  3.8.  Suppose  that  for  each  subset  y  of  [l,...,n},  the 
v 

mapping  L  is  norm-coercive.  Then  so  is  G. 

Proof.  Let  {xk}  be  a  sequence  with  [[xk[[  —  ®.  It  suffices  to  show  that 

for  any  two  partitioning  subsets  and  a of  {l,...,n}  and  for  any 

.  k, 

infinite  subsequence  (x  with 


[F(xki)  -  m(xki) 1 


k  k 

and  (BF(x  L)  +  (b  -  Cx  i))  =0 

*2 


for  each  k^,  the  subsequence  {|[G(x  ■*■) f[ }  diverges  to  ®.  For  any  such 
subsets  and  subsequences  we  have 

[G(xki)  =  l“2(xS]  and  G(x%  *  -<B  )’l[B_  ^  G(x%  -  (b  +  Ax%  ] 

1  2  2  2  2  1  1  2 


or  equivalently. 


/  k  \ 

k 

m 

GU  \ 

• 

/ 


<«„  o>  )_1  b« 

2  2  2 


/. 


B  B 

a  a  a  a 
2  1  2  2 


or,  k 
*  2(x  i)< 

k, 

^  2<*  i). 


Of  k  k 

As  \\t  2(x  *•)[[  —  therefore  so  does  [[g(x  i)||.  Hence  the  mapping  G  is  norm- 
coercive.  Q.K.D. 

Combining  Theorem  3.7  and  Proposition  3.8,  we  obtain  the  following 
result  concerning  the  (global)  existence  and  uniqueness  of  solution  to 
the  ICP(l). 


Theorem  3.9.  Let  m  be  convex.  Suppose  that 
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(14a)  For  each  subgradient  matrix  E  of  m,  the  matrix  given  by  (11)  is 
nonsingular  for  arbitrary  partitioning  subsets  a^,  a 2,  and  or3  of  Cl,...,n}; 

(14b)  For  each  subset  V  of  £l,...,n},  the  mapping  1?  is  norm-coercive. 

Then  the  mapping  G  is  a  global  homeomorphism  mapping  Rn  onto  Rn.  In 
particular,  the  ICP(l)  has  a  unique  solution  for  all  vectors  b  in  Rn. 

Proof.  Conditions  (14a)  and  (14b)  and  Theorem  3.7  imply  that  the  mapping 
G  is  a  local  homeomorphism  at  every  point  u  €  Rn.  The  first  conclusion 
now  follows  from  the  norm-coerciveness  theorem  in  Ortega  and  Rheinbolt 
[24,  p.  136-137].  For  the  second  conclusion,  it  suffices  to  combine 
Theorem  3.1  with  the  observation  that  the  conditions  (14a)  and  (14b)  are 
independent  of  the  vector  b.  Q.E.D. 

In  [22],  Megiddo  and  Kojima  have  derived  some  fairly  general  neces¬ 
sary  and  sufficient  conditions  as  well  as  several  more  specific  sufficient 
conditions  for  the  (global)  existence  and  uniqueness  of  solution  to  the 
NLCP(2).  As  pointed  out  in  the  Introduction,  if  the  matrix  A  is  non¬ 
singular,  the  ICP(l)  is . equivalent  to  an  NLCP.  This  observation  implies 
that  if  A  is  nonsingular,  then  the  results  of  Megiddo  and  Kojima  would 
provide  conditions  for  the  ICP(l)  to  have  a  unique  solution  for  all 
vectors  b.  This  nonsingularity  assumption  of  the  matrix  A,  although 
not  explicitly  made  in  Theorem  3.9,  is  nevertheless  a  consequence  of 
the  norm-coerciveness  of  the  mapping  -C,N(x)  where  N  -  {l,...,n}.  In 

v 

fact,  the  norm-coerciveness  of  the  mappings  L  (x)  is  equivalent  to  that 
of  the  principal  subfunctions  of  the  mapping  g  appearing  in  the 
corresponding  NLCP(2).  The  latter  is  one  of  the  conditions  needed  in 
Theorem  3.1  of  [22]. 
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4.  MONOTONE  ITERATES 


In  this  section,  we  study  the  ICP(A,b,m)  under  the  assumptions  that  m  is 
monotonically  nondecreasing  and  A  has  a  Z-splitting  (B,C)  with  C  nonnegative. 
Note  that  the  matrix  A  possessing  such  a  splitting  must  necessarily  be  a 
Z-matrix.  The  ICP(A,b,m)  with  A  and  m  having  these  properties  arises  for 
instance,  from  discretization  of  the  impulsional  control  and  optimal  stopping 
problems  mentioned  earlier  and  from  a  certain  generalized  LCP  (or  equivalently, 
the  problem  of  solving  a  system  of  piecewise  linear  equations)  solvable  as 
a  single  linear  program.  We  refer  to  Dolcetta  [13]  for  more  details  on  the 
former  application  and  briefly  discuss  the  latter.  Consider  the  generalized 


LCP 


(15)  M^x  +  q1  ^  0  for  i  -  l,...,k  and 


k 

TT  (M^x  +  q1)  =  0  for  j  -  1, 

i=l  3 


where  each  M*  is  n  by  n.  It  was  shown  in  Mangasarian  [21]  that  if  there 
1  k 

exist  Z-matrices  Y  ,...,Y  and  a  nonsingular  matrix  Q  such  that  the  two 


conditions  below  are  satisfied 


(16  a)  =  Y1 


k  .«  T  i  lk.  k  i 

{16b)  £  (s  )  Y  >0  for  some  (s  ,...,s  )  >  0  with  ^2^  s  >  0, 


then  a  solution  to  the  generalized  LCP(15),  if  it  exists,  can  be  obtained 
by  solving  the  linear  program 

T 

min  p  x  subject  to  M  x  +  q  >0  for  i  =  l,...,k 

k  A  T  4 

where  p  ■  2  (M  )  s  .  We  remark  that  the  result  stated  in  the  reference 


"T1  T 
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is  somewhat  more  general  than  the  one  presented  here.  (For  instance  Q 
needs  not  be  nonsingular  and  the  two  conditions  (16a)  and  (16b)  can  be  * 
relaxed.)  We  have  simplified  the  general  result  in  order  to  explain  our 
point  more  easily.  (See  Cottle  and  Pang  [11]  for  more  discussion.) 
Since-Q  is  nonsingular,- it  is  obvious- that  we  can  convert  the  generalized 
LCP  to  the  following  equivalent  one 

k 

Y^z  +  q1  >  0  for  i  =  l,...,k  and  'Tf (y^z  +  q1)  =  0  for  j  *  l,...,n 

i*l  J 

where  z  *  Q  ^x.  As  each  is  a  Z-matrix,  there  exists  positive  diago¬ 
nal  matrix  D*  such  that  >  Y^.  Hence,  with  A  =  Y^  and 

m(z)  =  max  [(I  -  (D^)  Y^)z  -  (D^)  qi]>  the  latter  generalized 
2  <  i  ^  k  x 

LCP  is  equivalent  to  the  ICP  (A,q  ,m)  where  A  is  a  Z-matrix  and  m  is  non¬ 
decreasing.  The  condition  (16b)  implies  that  the  feasible  set  of  the  ICP 
(A,q*,m)  is  bounded  below.  See  Cottle  and  Pang  [10]. 

Formally,  the  feasible  set  of  the  ICP  (A,b,m)  is 


X(A,b,m)  *{x€Rn:Ax+b>0  and  x  >  m(x)}. 


It  is  bounded  below  if  there  is  a  vector  s  €  R  such  that  x  >  s  for  all 
x  in  X(A,b,m).  The  following  result  shows  that  under  certain  conditions 
the  ICP  (A,b,m)  has  a  least-element  solution.  Its  proof  can  be  found 
in  Pang  [25],  ' 


Theorem  5.1.  Let  m  be  a  continuous  and  nondecreasing  mapping  from  R 
into  itself.  Let  A  be  a  Z-matrix.  Suppose  that  the  feasible  set 
X(A,b,m)  is  nonempty  and  bounded  below.  Then  X(A,b,m)  has  a  least- 


element  x*  which  is  the  unique  solution  of  the  program 


i 
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T 

min  ex  subject  to  x  €  X(A,b,m) 

where  e  is  any  positive  vector.  Moreover,  x*  solves  the  ICP  (A,b,m). 

A  mapping  f  from  Rn  into  RQ  is  inverse- isotone  if  f(x)  >  f(y) 
implies  x  >  y  for  each  x  and  y.  The  next  result  gives  a  sufficient  con¬ 
dition  for  the  feasible  setX(A,b,m)  to  be  bounded  below.  Its  proof  is 
easy. 

Proposition  5.2.  Suppose  that  for  some  subset  Y  of  {l,...,n},  the  mapping 
defined  in  Section  3  is  surjective  and  inverse  isotone.  Then  the 
feasible  setX(A,b,m)  is  bounded  below. 

Theorem  5.1  has  provided  a  least-element  interpretation  to 
Mangasarian's  result  of  solving  the  generalized  LCP  as  a  linear  program. 

For  more  discussion  on  related  results,  see  the  series  of  papers  [10,  11, 

17,  18,  19,  20,  26,  27]. 

If  (B,C)  is  a  Z-splitting  of  the  matrix  A,  then  the  mapping  F  is  not 
always  well-defined  for  two  reasons.  First,  F(u)  may  not  exist.  Second, 
F(u)  may  not  be  unique.  However,  if  the  set  X(u)  defined  in  Theorem  2.4 
is  nonempty,  it  has  a  least  element  by  Theorem  5.1  (see  also  Tamir  [32]). 

We  shall  let  F(u)  be  this  element,  so  that  it  is  well-defined  provided 
that  it  exists. 

Under  appropriate  conditions.  Theorem  5.1  shows  that  a  solution  to  the 
ICP(A,b,m)  can  be  obtained  by  solving  a  certain  nonlinear  program.  The 
next  theorem  describes  an  alternate  way  for  computing  a  (possible  different) 


solution  of  the  problem  under  the  same  set  of  conditions. 
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Theorem  5.3.  Let  n  be  a  continuous  and  nondecreasing  mapping  from  Rn  into 
itself.  Let  (B,C)  be  a  Z-splitting  of  the  matrix  A  with  C  nonnegative. 
Suppose  that  the  feasible  set  X(A,b,m)  is  nonempty  and  bounded  below.  Let 
u^  be  an  arbitrary  vector  in  X(A,b,m).  Then  with  u*^  as  the  initial  vector, 
the  sequence  of  iterates  [uk]  defined  by  :  uk+^  a  F(uk)  is  well-defined, 
nonincreasing  and  converges  to  a  solution  of  the  ICP(A,b,m). 


Proof.  By  induction,  suppose  that  {u^,u\  . . .  ,uk}  is  well-defined  and 

nonincreasing.  Moreover,  we  may  assume  that  each  u*  for  i  <  k  is 

k+1 

feasible  to  the  ICP  (A,b,m) .  We  now  consider  u  .  First  of  all, 

lc  tc  lc4”l  k*  ir|i 

u  6  X(u  ),  so  that  u  =  F(u  )  is  well-defined.  Since  u  is  the  least 

t  k  k+1 

element  of  X(u  ),  thus  u  >  u  .  Moreover,  we  have 


k+1  .  k.  ^  ,  k+1. 

u  >  m(u  ;  >  m(u  ) 

and 

_  k+1  ...  _  k+1  k+1  _  k.  ^  - 

k+1 

Hence  u  is  feasible  to  the  ICP  as  well.  Consequently,  the  entire 
sequence  [u  }  is  well-defined,  nonincreasing  as  well  as  feasible.  Since 
the  feasible  set  X(A,b,m)  is  bounded  below,  the  sequence  {uk}  converges 
to  a  vector  u*.  For  each  k  >  0,  we  have 

uk+1  >  m(uk),  Buk+*  +  (b  -  Cu^)  >  0  and  (uk+^  -  m(uk) )T(Buk+^  +  b  -  Cu^)  *  0. 


Hence  by  passing  the  limit  k  —  °®,  we  conclude  that  the  vector  u*  is  a  solution 
to  the  ICP  (A,b,m) .  Q.E.D. 
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This  is  the  first  of  two  papers  in  which  we  study  this  complementarity  problem 


via  an  implicitly  defined  mapping  F  which  depends  on  a  given  splitting  of  thi 
matrix  A.  In  the  present  paper,  we  derive  sufficient  conditions  for  the  problem 
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